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ABSTRACT 

This paper develops the mathematical relations 
for the '’Virtual Kinematic Model" as an improvement 
over the vectorial model developed earlier. The model 
solution described herein provides the most energy 
efficient means of cooling gears - that is, it requires 
the least pressure or pumping power to distribute the 
coolant on the tooth surface. Further, this nozzle 
orientation allows impingement to the root of the 
tooth if needed and provides the most cooling control 
when compared to into-mesh or out-of-mesh cooling. 

NOMENCLATURE 

a addendum of tooth 

d| impingement depth 

d x uncorrected impingement depth in origin plane 

ev of 6 = sec 4> = evolvent function of 0 

N number of teeth in gear 

n rotational speed 

P^ diametrical pitch 

P n normal diametrial pitch 

Pt transverse diametrial pitch 

R standard pitch radius 

Rb base radius of gear 

Rfcv base radius of virtual involute 

Rj radius at impingement point 

R or operational outside radius (for modified 
addendum) 

R 0 standard outside radius 


Rt transverse pitch radius 

R x uncorrected radius at impingement point (see 
Fig. 3) 

t arbitrary f 1 ight time 

tf jet stream flight time at impingement 
t w time to rotate through angle, 9 W 
Vg linear pitch line velocity of gear 

Vj jet velocity 

a jet i ncl i nation angl e 

Aa addendum modification 

d. P d 

S\ dimensionless impingement depth — 

cj v virtual roll angle at impingement point 

e 0 v virtual roll angle to outside diameter 

0 a vectorial angle of involute between point i 
and 0 (Fig. 2) 

9j involute of $\ - involute function at impinge 
ment point 

0iv virtual vectoral angle to impingement point 
9 0v virtual vectoral angle to outside diameter 
0 r angle across tooth space at base circle 
9 W total vectoral angle at impingement point 
Vj dimensionless jet velocity Vj/V g 

P b dimensionless base circle radius 

p. dimensionless impingement radius 
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p dimensionless outside radius 
r or 

P t dimensionless transverse pitch radius 

4> pressure angle at pitch radius 

$1 pressure angle at impingement point 


Rj after time (t = tf). From inspection of Fig. 1 , 
it can be seen that the radius of impingement is: 

Ri = R 0 - Vj tf (1) 

where Vj is the jet velocity. The time of flight of 
the oil jet through the tooth space is: 


1 v Pressure angle of virtual involute at impinge- 
ment point 

4 > n normal pressure angle 

<t>o pressure angle at outside radius 

<t>ov pressure angle of virtual involute at outside 
diameter 

$t transverse pressure angle 



The gear must rotate through the angle e w in 
the same time it takes the oil jet to travel through 
the tooth space to a depth of 6\ = Vjtf and impinge 
at R\. From Fig. 1 we see that e w = e 0 + © r + 9 j . 
The components of 0 W are developed below. 

The pressure angle, at the depth of impinge- 
ment, dj i s : 


Y helical angle 

Wg gear angular velocity 

Uj virtural jet angular velocity 


where 


cos 



and tan 4 >. 
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INTRODUCTION 


Rb is the gear base radius 


Cooling of lightweight high speed gears via the 
lubricating oil supply provided to the gearbox conti- 
nues to be one of the most difficult areas in the high 
performance mechanical power transmission design proc- 
ess. Since the cooling process is controlled by oil 
jet orientation and by the jet velocity relative to 
the pitch line velocity (PLV), these authors have con- 
sidered the development of a fundamental understanding 
of the oil jet primary impingement phenomena as essen- 
tial to understanding how the gear tooth cooling 
process works. In^ the early 1970’s when this work was 
started, a vectorial model" was used to describe the 
radial jet flow phenomena, (Akin, 1975), which com- 
pared reasonably well with the experimental data, 

(Akin, 1975), especially at pressures above approxi- 
mately 140 kPa (20 psi). Models were also developed 
for into mesh (Akin, 1983) and for out of mesh (Akin, 
1985) to determine the oil jet impingement depths for 
these conditions. The into mesh and out of mesh 
models indicate the oil jet impingement depths were 
limited and therefore maximum cooling could not be 
obtained using this approach. 

The 'Virtual Kinematic Model" was first mentioned 
by Townsend in 1980, but the full mathematical devel- 
opment of this model has not been published until now. 
This model has been generalized to include spur or 
helical gears with standard AGMA proportions as well 
as long or short addendums. Further, parametric limi- 
tations have been formulated to allow the gear analyst 
to know he is specifying his problem where a solution 
is possible. Also, the dimensionless parameter devel- 
opment has been provided to allow the gear engineer to 
examine a wide range of generalized solutions easily 
and to expose the limits mentioned above. 

KINEMATIC SOLUTION FOR MODIFIED SPUR GEARS 

The jet stream flight time starts (t * 0 ) at 
position " 1 " at the intersection of the top land and 
trailing profile of the leading tooth, as shown in 
Fig. 1 . Then the terminal end of the jet passes 
through the tooth space where it impinges at " 1 " on 
the leading profile of the trailing tooth at radius 


Rj = Rq - d j 

so that 0 j = inv $j = tan <$>j - 

Also, the angle across the tooth space at the gear 
base circle is: 


r = N 


2 i nv 


2 ( tan $ - <j>) 


(4) 


where <t> = cos 1 ^ R 5 /R) is the pressure angle at the 
standard pitch circle R. Further, the pressure angle 


at the outside radius R 0 is: 


f%\ 

R 2 

and tan * = 

IV 0 J 

■V i 
[ R b J 


1/2 


(5) 


so that 


and thus 


e 0 = inv *0 = tan *0 


e w = e 0 + e r + 9 i as 5tated earlier, so that the time 
it takes the gear to rotate through 6 ^ is 

e w 30 e w 

*w * Z — sec; < 6 > 

g 

where n = rev/min. 

Equating the times of flight t. and rotation t ■ 

i w' 


d. 30 0 
i w 


V. 

J 


irn 


or 


R o- R i 

V, 


(7) 


v j u g 

If d 1 is given, then the required V. can be found 
expl i ci t ly as 


2 



V 


j = 



ird . n 


From here we can develop the remainder of the virtual 
(8) involute functions including 9 0V and 4> ov as fol- 
lows (see Fig. 2): 


note, that d. is a function of <$>.(d, = f (<M) 
l l i n T i 

If \/^ is given and the resulting d. or R. 
is desired, we must iterate until R. , as fn($. , V.) , 

is satisfied. In dimensionless form, where: 

V. = V./V g and 6, = di P d /2 

d, 45. 

v j - R o e w = [<N ♦ 2)6 W ] (9) 

GEARS WITH MODIFIED ADDENDUMS AND AN OIL JET DIRECTED 
AT AN ARBITRARY INCLINATION ANGLE a FROM THE GEAR 
CENTER 

A more general solution can be provided to this 
problem if we assume that the jet is pointed away from 
the line through center of the gear by an arbitrary 
angle "a" as shown in Fig. 2. This angle is defined 
to be positive in magnitude when pointed in the direc- 
tion of wg as shown in Fig. 2. Modified addendum 
lengths or center distances are accommodated by the 
use of the parameter AR 0 = ANg/2P n = Aa, so that 
R or - R o + Aa. Obviously when Aa =0, R 0 = R or . In 
this paper we will generalize and use Aa, R 0 , and 
R or in al 1 equations . 

For helical gears we let: * N ( 2P n cos y) , so 

that, Rq — R{- + a and Rq^* — R^ + a + Aa, where 
p t = p n cos Y and the subscript t refers to the 
transverse plane (normal to gear axis.) This problem 
has been formulated by assuming that the jet terminal 
end (or head) follows a trajectory in mesh with a 
"virtual involute" as shown in Fig. 2. The virtual 
involute is: defined by the radius of the base cir- 

cle Rbv. where 


R bv = R or s1n «• 


V. 

J 


. R u , 
J bv’ 


or 


“j ■ R 


V. 


or 


si n a 


( 10 ) 

As can be seen by inspection of Fig. 2, "u." is 
the theoretical angular velocity of the virtualjet 
involute that produces the proper locus for the term- 
inal end of the jet stream line at the actual jet 
velocity Vj, subsequent to passing the point 1 at 
t = 0. The virtual involute is used to develop the 
mathematical relationships needed to formulate the 
time of flight to impingement of the jet stream on the 
trailing tooth profile at when the gear angular 

velocity is wg. Note that in general cjg * w. since 

the virtual base circle is not fixed on the gear blank 
(or wheel ) . 

The writer has selected the virtual involute gen- 
eralized roll angle e v as the parameter to interre- 
late the virtual involute for the jet stream with the 
rotating gear geometry. Thus the virtual roll angle 
e ov at the gear outside diameter can be calculated 
from the fundamental relationship (see Fig. 2): 


R bv e ov 


R 0 r cos a * so thdt c 0 v 


R or COS a 
R bv 


R or cos * _ _L_ 

R_ sin a tan a 
or 


cot a <0 < a < j - ) 


( 11 ) 


7 1 17 7 \ 17 

e ov = ^ ev ^ - 1' or ev e = (e^ v + 1) ' where by 

definition (Vogel, 1945): 


ev 6 = W 50 that E ov - [(R or /R bv )2 - 1]W2 

avoids the use of a functions. The virtual 
vertorial angle or polar angle 9 qv at the O.D. is: 


9 ov = [R or /R bv )2 - ]]UZ - tan ~ 1t(R or /R bv >2 ' n’'* 
(rad.) and the virtual pressure angle is 

obtained from: 


i s inv 9 
r ov ov 


sec' 1 (ev 9^ ) 
ov 


cos 


-1 


(R, /R ). 
bv or 


( 12 ) 


From the above we get e =9 + <±» as a check. 

J ov ov Y ov 

Also, a = a = -rr/2 - 4> (rad) as a check) 
ov T ov 


ov 


R^ COS a 
or 

R bv 


= cot a 


The virtual roll angle e. 


1 v 


diameter is developed as follows: 


(0 < a < |-> < 13) 

at the impingement 


c . — c -e. = co • t ^ = 

0-1 ov iv j f 




si n a 


(14) 


so that e.(t^) = e QV - ojj t ^ = (R Qr cos a - V^t^)/ 
(R Qr sin a) = cot a - Vj t ^/ ( R Qr sin a). 

Note from Fig. 2 that V^t^ = d./cos a, thus when V^ 
and/or t^ are unknown: 


e 1v (d l> ■ 


R or cos 


d i 

COS a 


d s 


R A „ sin a 
or 


= cot a - 


R sin a COS a 
or 


and e . = e - e. 
o-i ov iv 


d: 


2d, 


COt a - COt a 


R Qr sin a cos a ~ R Qr sin 2 a 


(15) 


Again the other virtual involute functions at the 
impingement diameter are; starting with the vectoral 
angle : 


i v 


(R - d,) 2 

1/2 

(R - d, > 2 

or i 

- tan -1 

or 1 

2 ' 
R bv J 

2 " ‘ 
R bv J 


(rad. ) 


and the virtual pressure angle 4 >- v at R. : 


( 16 ) 
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? iv 


= cos 


-1 


bv 

R or " d ij 


( 17 ) 


Further noting that 


iv 


(R or-V‘ 


bv 


1/2 


r iv * v as a c ^ck, an ^ 


(18) 


a iv = r- ^v (rad - } and R i - R bv (E iv + ° 


However, we cannot use dj to calculate e a 
except as an approximation because the angular origin 
for 6 a is along the radius vector R x (Fig. 3) pass- 
ing through the time origin at "O'* while dj is meas- 
ured along the Rj radius vector terminating at point 
"i” on the leading profile of the trailing tooth, 
which is the terminating vector position for the’ 
angle 0 a ; thus (see Fig. 3): 

R x » Rj cos © a = (R or - dj)cos 0 a (22) 

dx = Ror ~ Rx = Ror - (Ror - dj)cos 9 a (See Fig. 3) 
Noting that: 


as a check. Thus the time of flight can be checked 
f rom : 


d x tan a = R x tan © a 
it can be shown that: 


( 23 ) 


t - _2nl 
r f ‘ " R 


d; 


Ij R or sin a cos a[Vj/(R or sin a >] 


d . 


Vj cos 


(sec) (19) 


where 


dj = (R 0 + Aa) - Rj * R or - Rj 

Theoretically, the virtual involute is a dummy 
device for mathematically describing the time motion 
of the terminal end of the jet stream line. As such, 
the jet stream line is the line of action between the 
virtual involute and the oil jet. Thus the coincident 
points "i" at radius rj on the virtual involute 
(v.i.) and Rj on the gear tooth (g.t.) define the 
point of oi 1 impingement at time tf (v.i.) = t w 
(g.t.). This coincident instant in time is the result 
of the simultaneous motions of the (v.i.) and (g.t.) 
from time t = 0 when the (v.i.) is at its outer 
position 2 at the virtual base circle (Fig. 2) while 
the trailing edge of the leading gear tooth is at 
position M 0" where the jet line crosses the 0.0. and 
the (v.i.) is such that the rotation of the (v.i.) at 
wg without the relative rotation of (w. - cog) through 
time t w would place the (v.i.) origin^at position 3 
at the virtual base circle. Thus the (v.i.) is pre- 
ceded to rotate faster than u. by the amount 
(w j - “g> such that it rotates at w g + («. - Wg ) = «. 
in keeping with the jet velocity Vj whicfi places the 
(v.i.) at position 4 at the virtuarbase circle and at 
i on the leading profile of the trailing tooth as 
shown in Fig. 2 

It is now necessary to involve the gear tooth 
geometry. The solution will involve the assumption 

Tj = Rj = R or - dj 


where 


dj = V.t f cos a (20) 

Thus from the geometry of Fig. 2 it can be seen that 
©a can be approximated from: 


tan 9 a ■ r 


Vjt f sin a 


d. tan a 


or 


Vjt f COS a 


Or 




( R Qr - Rj >tan a 

R. 


( 21 ) 


dx * Ror cos a - [(Ror cos a)^ - 

dj(2R or - dj)] 1/2 COS a. (24) 

so that, 6 a = tan- 1 [d x tan <x/(R or - d x )]. 

Then as a check on the calculated value of d x 
above: 


R - d 

d r . -o r ... _ x 
i or cos 9. 


which should be the same as before. 


(25) 


If it is noted that R. = R Qr - d. then 
(see Fig. 2). 



( 26 ) 

Note that <(>. is a function of Ri and dj making an 
explicit solution for (Rj or dj) Impossible when Va 
is specified leaving dj or Rj as the dependent J 
variable. Also. 
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r 


TT_ 

N " 


4Aa tan 4> 
_n 

R 


29 n 


(27) 


where: = inv 4> n = tan 4 > n - ^ for spur gears and 

e r = ir/n - 2Aa tan <t> t /R t - 20 t for helical gears 


where 


©j. - inv = tan <|> t - and tan = tan <t> n sec y 



Now the gear angl e of rotation, ©w, between points 1 
and i in Fig. 2, can be calculated from 


© = © 
w o 


+ e i 


(30) 
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so that 


( 38 ) 


e w 

t = — 

W 0 >_ 


30<e o ♦ e r 


e, ♦ e a ) 


where 


irn 


n = rev. /min 


sec 


( 31 ) 


It Is now possible to equate the time of flight 
tf of the oil jet to the time of gear rotation t w , so 
that 


where 


t 3 

W Ug 


9 w = 6 o + 6 r + 6 1 + 8 a 


so that angular geometry expressed in fundamental 
terms of N, AN, <f>, a and 6. or are; 

starting with the vertorial angle 0 (see Fig. 1 

a 


yg-t-) 


3oe., 


rrn 


(V. COS a) = t f <V -'- ) 


( 32 ) 


Since 9. is fn(d.), the impingement depth 

cannot be solved explicitly but must be iterated 
numerically. Thus, we solve for so that an 

explicit solution is possible as follows: 


ird . n 

w = ! 

J 300 COS a 


0 t COS a 
w 


( 33 ) 


It would be desirable to graph the results of 
these formulas as in Figs. 4 and 5. Thus the 
independent variables will be normalized to provide 
dimensionless solutions to the above results. The 
principle parameters used are Vj and 6\ so that 
we will express them In terms of the dimensionless 
parameters , S\ and N as: 


v 3 . 

3 V g ‘ “- R 


1 


g or 


8 w R or cos “ 


(34) 


and 


d, d, d. P 


"1 "i n , . 2S i 

’J ’ v j“ , g r 'or°l = H.D. = TTT = ~~r and d i ■ — 

n n 


V . - v.cdR S. = 


( 35 ) 


where W.D. is the whole depth. 

Thus the time of flight (t^) from the trailing 
edge of the leading tooth top land at radius R Qr 
to the radius of impingement R. on the leading 
profile of the trailing tooth is: 


2S i 


f ” v.u PR COS a 
j g n or 


(36) 


and in terms of N, using R = (N + 2 cos ^ + 2An)/ 

(2P cos y): 
n 


46 ^ cos y\t 


(37) 


*f = [v^cOg(N + 2 cos y + 2AN)cos a] 

Then the time of rotation t from the trailing 

w 

edge of the leading tooth top land at radius R Qr to 


to the jet line when it intercepts the impingement 
point at radius R 

trai 1 i ng tooth 1 s : 


point at radius R Qr on the leading profile of the 


tan 0^ = d. tan , so that 


0 - tan" 1 

a 


46. 

N + 2 

cos y 

Next the impingement 

angle 

0. = tan <$>. 

, - A. 
i Y i 

.(4 



\ p b 


J/2 


- cos 3 

lP i. 


where 


thus 


and 


46. cos t|i 

p i = 1 N + 2AN + 2 COS «|> 
N COS 

p b N + 2AN + 2 cos y 


((N + 2 cos y + 2 AN) - 46 . COS i|r) 


N cos 4> t 


p = 1 
F or 


Now this is valid only when (N is given): 


s i i 


N + 2AN + 2 COS y - N cos 4>. 


4 COS y\i 


or when S\ is given 


N > 2 


( 26 . - 1 ) COS I|r - 


AN 


1 - cos 4> f 


6. (max) 


= N(min) 


the tooth space angle 0 r at the base diameter is 
it 4 AN 


e r ■ r- T tan * t - 29 t 


where 


0j_ = tan - 4> t - 


^ - 1 


1/2 


-1 3 
cos — 


(39) 


i s : 


(40) 


(41 ) 


(42) 


( 43 ) 


( 44 ) 


( 45 ) 


(46) 
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where 


Pf = 


N 


and 


t “ N + 2 cos 4 / + 2AN 


!t _j 

P b " cos <t> t ■ sec ♦t 


and finally the involute angle 8 0 to the outside 
diameter is: 


i nv * = tan d> -6 

or Y or y or 


1/2 


cos“ (p b ) 


(47) 


again it is possible to equate the times of flight t^ 
and rotation t w , now in terms of dimensionless param- 
eters , so that 

e. 


4 5. COS y\t 


Vj ( N + 2 COS \\> + 2AN)cos a 


(48) 


Therefore, if the dimensionless depth s. < $. 

(max) is given the dimensionless velocity is cal- 
culated from: 


V j " 


4$. COS 4 / 


0 W <N + 2 cos 4 / + 2AN)cos a 


(49) 


jet nozzle. Also the effect of jet velocity on depth 
S\ is shewn in Fig. 5. The curve presented here is 
for a specific gear design: pressure angle 4 = 20 °, 

number of teeth N = 28, inclination angle a = 47° 
and dimensionless jet velocity Vj/V g = v . = 0.6468. 

Figure 6 (Akin, 1975) shows the theoretical 
results compared with the experimental data gathered 
in 1974. The experimental data did lie substantially 
below the prediction from the mathematical models at 
low oil pressures. It is assumed, as also mentioned 
earlier, that this is the result of ( 1 ) the windage 
effect on a fan-jet at low jet velocity and ( 2 ) pres- 
sure loss in the nozzle. As the jet velocity is 
increased the experimental data approaches the curves. 
The windage study conducted in Akin (1975) involved 
predicting the trajectory of the oil droplets after 
they passed into the tooth space between the gear 
teeth. Under these circumstances the windage affect 
is considered negligible except for small droplet 
sizes around 0.001 cm (0.0004 in.). 

CONCLUSION 

This paper develops the mathematical relations 
for the "Virtual Kinematic Model" as an improvement 
over the vectorial model developed earlier. The model 
solution described herein provides the most energy 
efficient means of cooling gears-that is, it requires 
the least pressure or pumping power to distribute the 
coolant on the tooth surface. Further this nozzle 
orientation allows impingement to the root if needed 
and provides the most control when compared to the 
into-mesh or out-of-mesh cooling. 
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FIGURE 3. - ENLARGED DETAIL OF KINEMATIC MODEL TO SHOW GEOME 
TRY CORRECTION FOR di AS A RESULT OF 0„. 
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DIMENSIONLESS JET VELOCITY \)j 

FIGURE 5. - DIMENSIONLESS 6j FOR A GIVEN DIMENSIONLESS 
JET VELOCITY Vj (AT OPTIMUM VALUE FOR ALPHA) a = 47°, 
PRESSURE ANGLE, 20°; NUMBER OF TEETH. 28; 8 DIAMETRIAL 
PITCH. 
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FIGURE 6. - CALCULATED AND EXPERIMENTAL IMPINGEMENT 
DEPTH VERSUS OIL JET PRESSURE AT 4920 AND 2560 RPM. 
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